EXISTENCE AND UNIQUENESS OF AN INVARIANT 
MEASURE FOR A CHAIN OF OSCILLATORS IN 
CONTACT WITH TWO HEAT BATHS 



PHILIPPE CARMONA 

Abstract. In this note we consider a chain of N oscillators, whose 
ends are in contact with two heat baths at different temperatures. Our 
main result is the exponential convergence to the unique invariant prob- 
ability measure (the stationary state). We use the Lyapunov's function 
technique of Rey-Bellet and coauthors [TT1 [H El [HI 111 IS] , with different 
model of heat baths, and adapt these techniques to two new case recently 
considered in the literature by Bernardin and Olla [2] and Lefevere and 
Schenkel 0. 



1. Introduction 

We consider a chain of coupled oscillators whose dynamic is described by 
the Hamiltonian 

H = E \pf + V ^)+ E U(q i+1 - gi ), 
l<i<N l<i<JV-l 
where qi and pi are the position and the momentum of the oscillators. We 
assume that the potentials U, V have the following properties: 
(HI) : Growth at infinity U and V are C°° and there exists real constants 
I > 1, k > 2, k > I, ak > 0, bi > such that: 

lim X~ k U(Xx) = ak\x\ k , lim \ 1 ~ h U'(\x) = kdk\x\ k ~ l sign(x) 

A— >+oo A^+oo 

lim \~ l V{Xx) = b l \x\ l , lim X l - l V {\x) = l ai \x\ l ~ l sign{x) 

A-— >+oo A^+oo 



(H2) : Non degeneracy The interaction potential is non degenerate : for 
any q € R there exists m = m(q) > 2 such that d m U{q) ^ 0. 

We want to stress the fact that k > I, that is near infinity the interaction 
potential U dominates the pinning potential V. Typical examples of such 
potentials are: 

• The harmonic chain : U(x) = \x 2 , V = ax 2 , with a > 0. 

• The Fermi-Pasta-Ulam chain : U{x) = \x 2 + V(x) = ax 4 , with 
a > 0. 
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The harmonic coupling with weak pinning : U(x) = \x l 



V{x) 



(1 +x 2 )2. 

The two ends of the chain are in contact with heat baths at temperatures 
T\ = Ti > T r = Tjy. The interaction between the heat bath and the particle 
is modeled by a Langevin process (an Ornstein Uhlenbeck process) at the 
corresponding temperature. Therefore the dynamics are described by the 
following system of stochastic differential equations: 



(S) 



1 



( d qi (t) =Pi{t)dt 
d Pi (t) = (-d q .H 

[d Pi (t) = (-d qi H)dt 



Pi ) dt + y/TidBiit) 



(i = l,N) 
(2 < i < N — 1) . 



where B\,Bn denote two independent standard Brownian motions. 
In the unpinned case, that is V = 0, the dynamics are translation invariant 
(ie under the action qi — > qi + C) , and there is no hope of finding an invariant 
probability measure. It is thus natural to consider the interdistances Ti = 
Qi+i — <7i, for 1 < i < iV — 1, and the system 



(S>) 



Pi) dt 
1 



dn{t) = {p i+ i - 
d Pl (t) = (U'(n) 

dp N (t) 

{ dp i (t) = (U'(r i )-U'(r i _ 1 ))dt 



pi)dt + y/TidBi(t) 
1 



(-WfXN-i) - -Pn) dt + ^/T N dB N (t) 



(2<i<N—l) 



Observe that the state space of the stochastic system is x 
pinned case, and x = ^ 2N ~ 1 in the unpinned case. 



R 2N in the 



To state the main result, we need to introduce the Lyapunov function W 



8H 
space 



e"" , for a 9 > to be chosen later, and the corresponding weighted Banach 



n = {f-. x 



continuous, 



< +00} , 



with 
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sup^M 

s 7 w(x) 



Theorem 1.1. There exists 8 > such that: 

(i) For every starting point, the system (S) (resp. {S')) has a unique 
solution, defined for all times t > 0. 

(ii) The corresponding semigroup (Tt)t>o has a smooth density. 

(Hi) The Markov system has a unique invariant probability measure tt, which 

is absolutely continuous with a smooth density, 
(iv) The semigroup converges exponentially fast to the invariant measure 

(and is therefore ergodic). More precisely, there exists a constant C > 

such that 

\T t f(x) - tt(/)| < C\\f\\ w e- ct W(x) (Vx 6 Xl /eH) 

The proofs we give are similar to the original ones, as they are stated in the 
survey paper [10J, or in the original papers by Rey-Bellet and Thomas [8| [T3] . 
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and Eckmann, Pillet and Rey-Bellet [HE]. However we believe that this note 
is an interesting contribution to the subject, because of the following: 

• The model of heat baths is the simplest possible, an Ornstein-Uhlenbeck 
process. We do not introduce the auxiliary variables of [SI EH HI El E] • 
However, we do not provide a nice physical interpretation of such a 
model of heat bath. 

• The assumptions on the potentials are slightly relaxed : we do not 
require that the pinning potential to be at least quadratic, that is 
I > 2. Observe that we do not reach the level of generality of Eck- 
mann and Hairer [3] , where only growth bounds on U, V and their 
derivatives are required. 

• In [HI [131 HI El E] , the authors use a scaling argument based on two 
ingredients : 

— Continuity of solutions of ordinary differential equations with 
respect to initial conditions and coefficients. 

— A tracking lemma which says that at high energy the system 
is nearly deterministic (one can almost forget the heat bath 
influence). 

We replaced both arguments by a single one : continuity of sto- 
chastic differential equations with respect to initial conditions and 
coefficients. This is not only a technical shortcut, but also enables 
the study of models of chain where in addition to heat baths, you 
have a specific noise that you need to keep even at high energy since 
it is a feature of the model. 

• The first of these models was introduced by by Bernardin and Olla [2] . 
There is, in addition to heat baths, a random exchange of momen- 
tum between neighbouring atoms. The model at infinite energy is 
not deterministic but stochastic (see the system |E°°1 on page [T3j) . 

• The second model we study is the chain introduced by Lefevere and 
Schenkel [7], where the noise is highly non standard : in Fourier 
coordinates the momenta are in contact with heat baths at the same 
temperature T and positions are also coupled to the heat baths. 
It is not a straightforward application of results of Rey-Bellet and 
coauthors since we need to show that the energy dissipation on the 
momenta is the leading term in the upper bound of (|8.ip . Alas 
we are not able to generalize the results to non quadratic pinning 
potentials V. In particular we do not know how assumptions of the 
type C\ < V" < C2 could be translated in Fourier coordinates. 



with locally Lipschitz coefficients b : R n — > R n , a : M. n —> M nxmi B a 
Brownian motion in K m , we can only be sure, given a starting point x, of the 
existence of a local solution, defined up to an explosion random time C(w). 
However, if there exists a Lyapunov function W, ie a function such that 
W > 1 and the level sets {W < A} are compacts, that satisfies LW < CW 



2. Non explosion 
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for a constant C > 0, with L the formal generator of the diffusion X, that 
is the second order differential operator 



CCr*l 



then (see Theorem 5.9 of [9]) the lifetime is £ = +oo, that is X is defined 
for all times t > 0, with semigroup (Ti)t>o such that : 

This ensures that (Tt)t>o can be extended to a strongly continuous semi- 
group on H. 

For the pinned chain of oscillators, we have 

L = A + £ R , with A = Y J d Pl Hd qi -d qi Hd Pt 

i 

the Hamiltonian generator and Cr the heat baths (reservoirs) generator 

C R = \ E (-Pi^+T^). 

i=l,JV 

The energy H is conserved by the Hamiltonian dynamics : AH = 0. Since 
A is a first order differential operator 

AW = e m 6A{H) = . 

We shall make the assumption 

0<9< -J—— 

Then, 

C R H = ^(T l -p 2 1 +T N -p 2 N ) 

LW = C R W = ±6W((Ti + T N ) - {p\+p 2 N ) + 9{T 1P \ + T N p%)) < ^(T x + Tjy) TV 

For the unpinned chain, we obtain exactly the same upper bound, since we 

still have AH = AW = 0. 

Hence, part (i) of Theorem II .11 is established. 

Observe that since V(x) and U(x) go to infinity as |x| — > +oo, H is bounded 
below, therefore we need to multiply e SH by a constant C in order to have 
W = Ce 8H > f. We shall forget about this constant in the sequel since it 
will not change any proof. 

3. Hypoellipticity 
The generator of the pinned system can be written 

L = Xq + X\ + Xjf , 
where Xi are first order differential operators 

X = A - -{pid pi + p N d PN ) , Xi = \j^-d Pi (i = l,N), 
and A = Y^i 9 Pi Hd qi — d qi Hd Pi is the Hamiltonian generator. 
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Let £ be the Lie algebra generated by the vector fields 

{ X i}i>l> X j]}o<i,j>{[Xi> [ X 3> X k]]} <i t j tk 

Assume Hormander's hypoellipticity condition, that is that £ has full rank 
at every point x. Then (see e.g. Corollary 7.2 of [§] and the references 
therein, or Hormander's original paper [6]), 

• The semigroup Tj has a smooth density, i.e. 

T tfix) = j p t (x,y)f{y)dy 

with (t,x,y) — ► Pt(x,y) smooth. 

• The semigroup Tt is Strong Feller : for t > 0, Tj sends measurable 
bounded functions into continuous functions. 

• The invariant measures, if they exist, also have a smooth density. 

Hence, we only need to verify that £ has full rank, in both pinned and 
unpinned cases. 

Pinned case £ contains Xi,Xn so it contains d pi and d Pn . But, 

[d Pl ,X ] = --d Pl + [d Pl ,A] = --d Pl +d qi 
Therefore, d qi G £. Since, 

[d qi ,X ] = [d qi ,A] = - d qiQiH d Pi = -d qiqi Hd pl - d qiq2 H d P2 

i 

we have d qiq2 H d P2 £ £. We iterate this procedure, by computing [d qi , d qiq2 H d i 



P21 

until we obtain that d q ^i q2 Hd P2 = (-l) m ~ l d m U{q 2 - qi) d P2 G £. By the 

non degeneracy assumption, this implies d P2 G £ and [d P2 , Xq] = q 2 G £. By 
induction, we establish that for all i, d Pi and d qi are in £. 

Unpinned case Only the Hamiltonian generator is different : 

A= Yl (Pi+i~Pi) d n+ (^(^)-^(^-i))^ l +^ / (n)pi-C/ , (r iV -i)a pjv . 
l<i<JV-l 2<i<AT-l 

By assumption d Pl G £ and 

[d Pl ,X ] = --zdpi ~ ^i 

Hence d ri G £ and 

[d ri ,x ] = -U"( n )d V2 + U"( n )d pi 

Consequently, U"(ri)d P2 G £ and 

[d n ,U"( ri )d P2 ] = d 3 r3 U( ri )d P2 G £ 

and we obtain by induction that d^V {r\)d P2 G £. Therefore, d P2 G £ and 
since 

[dp2i X o\ = 9 ri d r2 

we have d T2 G £ and by induction we show that £ has full rank IN — 1. 

Remark 3.1. Observe that we only need one heat bath to ensure hypoellip- 
ticity. 
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4. The Control Problem 

The way chosen to establish uniqueness of the invariant measure, if it ex- 
ists, is to prove that the semigroup (T t )t>o is irreducible, that is for every 
x, every t > 0, every non empty open set A: T t (x,A) > 0. The most fre- 
quently used tool to establish irreducibility is Stroock-Varadhan's support 
theorem|14|: let X be the solution of the stochastic differential equation, in 
the Stratonovitch sense, 

8X t = b{X t )d t + a(X t )dB t 

Let U = {u : R + -> R m : u(0) = 0,u' G £ 2 (0,t)Vt} be the set of controls 
lying in the Cameron-Martin space. Let (Cf, iCTVU ) be the controlled ordinary 
differential equation 

(4.1) ^ = b{x) + a(x)u{x) 

Let x(xQ,u,t) be the maximal integral curve of the controlled equation, 
passing through xq at time t = 0. Let A(xo,t) (resp A*(xo,t) ) be the 
accessibility set (resp the strong accessibility set) at time t > starting 
from xq- 

A(x, t) = {x(xq, u, s) : u G U, < s < t} , A*(x, t) = {x(xo, u,t) : u € 

The the support theorem states that the support of the measure T t (xo, .) is 
cl(A*(xQ,t)) the closure of the strong accessibility set. 

In the general case, establishing that cl(A*(x ,t)) = x(= ^- 2N or M 27V_1 ) 
is not a trivial thing, and we refer to Eckmann,Pillet and Rey-Bellet[5], 
Theorem 3.2 for a proof. 

However, for the harmonic chain, U(x) = ^x 2 and V(x) = \ax 2 , we can 
provide a simple proof, that we learnt from Frangois Laudenbach. Indeed, 
the control system is linear. First since the matrix o{x) is constant, there is 
no difference between ltd and Stratonovitch integrals. Second, let us write 
the control system, for the pinned case: we incorporate the constants \/J{ 
into the two controls m,ujsr so that we obtain 



(C) 



( d qi (t) = Pi (t)dt 
d Pi {t) = (-d Qi H - - Pi ) dt + Ui (i = 1, N) 



\ d Pi {t) = i-d qi H) dt (2 < i < JV— 1). 

This can be written, if x = (p, q) G M. 2N as 

dx . „ , . 
— = Ax + Bu(t) 

with A a n x n matrix and B a n x 2 matrix (n = 2iV). Kalman's criterion 
states that for every t > and x G x = we have ^4*(rco, t) = M n as soon 
as the smallest vector space S containing the image of B, and stable by A, 
is the whole space \ = ^ n ( see e -g- Wonham's|15|). 

Here the image of B is spanned by {ejv+i, &2n} where the are the canonical 
base of M. n . Assume for example that V(q) = \oiq 2 and U{q) = \q 2 ■ then 
the control system is 
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(C) 



dqi 
dt 

dpi 
dt 
dpN 

dt 

dpi 
dt 



Pi{t) 

q 2 ~ qi+ aqi + ui 

~(qN - qN-i) + aq N + u 2 

{qi+i + qi-i - 2%) + aqi 



(2< i< N-l). 



The matrix A can be written as blocks 

A = 



A 

1 

Therefore we only need to check that the smallest vector space containing 
ei and ejy, and stable by A is R N itself. This is obvious since we have 



A = 



f-l + a 1 
1 -2 + a 1 











\ 



-1 + a) 



5. Scaling 



We shall scale both positions qi and moments pi and use continuity with 
respect to parameters, for solution of stochastic differential equations, to 
prove that the system is attracted to the compact sets where the energy 
H stays bounded. More precisely, given X(t) = (q(t),p(t)) solution of (S) 
starting from x = (g(0),p(0)), and E > 0, we shall consider 



X E {t) 



' pf{t)=E-^p i {E 1 k-h) 



qf{t)=E~t qi (Et-H) 



Of course = pf and for 2 < i < N — 1, we have 



dt 



dtf 

dt 



-d qi H E (q E (t)) 



with the Hamiltonian 

H E (p, q)=J2 hi + l V ^qi) + E ^U{E^{q t+1 - • 



l<j<Af-l 



For i = 1, iV we obtain 



dp 



(t) = E-*p(0) ~ J Q (E^-pf(s) + d qi H(q E (s))) ds + E^^Bf{t) 

where B E (t) = E*~2k B{Ek~n) is a standard Brownian motion. In a 
nutshell X E is the solution, starting from x E = (E~ip{Q), E~ *q(0)) of the 
stochastic differential system 
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(S E ) 

( dqi(t) = Pi (t) dt 

d Pi (t) = {-d qi H E - -JLlpJdt + En-iy/TidBiit) (i = 1, AO 

, dpiit) = {-d qi H E ) dt (2 < i < N - 1) . 

Observe that the scaling is such that if H(x) = E then He(x e ) = 1. 

The main ingredient of the proof is the convergence of a solution of (S E ) to 

a solution of the limit system, where the noise has disappeared, 



(S°°) 



with 



( dqi(t) =Pi(t)dt 

dPi(t) = {-d^H^ - i l( fc= 2)Pi) dt (i = 1, N) 

{dp i {t) = {-d qi H 00 )dt (2<i<N-l). 



Hoo(p, q) = ^2 ^ + 1 (k=i) b k\qi\ k + a k\n+i - n\ k - 

i Ki<N—l 



We shall establish now that the system (S°°) is non degenerate, that is that 
the chain remains still when both ends are still. 

Lemma 5.1. Assume I = k. Let (q(t),p(t)) be a solution of (S°°), starting 
from x such that H QO (x) = 1. Then, for any r > 0, 



[ (p?(s)+p^(s))ds>0. 

J 



Proof. Assume that Jq(pKs) +p 2 N (s))ds = ; we are going to obtain a 
contradiction. Let / = [0,r]. Then, by assumption, p\ = on / so = 
pi = and qi = c\ is constant on /. We have, on / 

= -Jr = -( Hk=2)\p\ + dq.Hvo) 

therefore 

= 6 fc |ci| fc ~ 1 sign(ci) + a k \q 2 - ci| fc_1 sign(g 2 - ci) 

hence q 2 = c 2 is constant on /, and by induction we obtain that pi = and 
qi = Ci are constant on /. If we write H OQ (p, q) = Y^i \v\ + g{q) we obtain 
that c = (ci, . . . ,cjv) is a solution of Vg(c) = 0. Since 5 is convex (k is 
even), c is a global minimum of g, which is absurd because g(0) = and 
g{c) = #00(0, c) = 1 by assumption. □ 
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Remark 5.1. For the unpinned case we obtain exactly the same result, for 
the limit system 



dri(t) = {pi + i - Pi) dt 



ijfc-l 



dpi(t) = {ka k \n 
dpN(t) = (-ka k \r N -i 



1 



^ l(k=2)Pl)dt 



2 

fc-i 



dpi(t) = kak(\n 



2 

lfc-1 



l(fe=2)PAr)rf* 
) (it (2 < i < iV - 1) . 



If H oa (x) = 1 and X a solution of (S'°°), starting from x, then /J"(Pi(s) + 
p^(s))cis > 0. 

From the preceding Lemma, we shall derive the asymptotic result: 

Lemma 5.2. Assume that (X n (t))t>o is a solution of ($) starting from x n , 
with E n = H{x n ) — ► +oo. Then there exists a subsequence, {x nh ) k such that 
for any C > 0, to > 0, 



lim P Xn k 
k— >+oo 



exp 



to 



C (pi(s)+p 2 N (s))ds 



0. 



Proof. The process (X,f n (i)) t >o is a solution of (S En ) starting from x^ 
with H En (x^) = ^-H(x n ) = 1. We build, on the same filtered probability 
space, that is with the same Brownian B, a family of processes (X^ n \t))t>o, 
solution of (S En ) starting from x^ n \ 

Assume first that k = I. Since He h (x^) = 1, the sequence x n remains 
in a compact set and we can extract a converging subsequence, which we 
shall still denote by x^ n \ Hence we have 

x (n) _^ x (oo) 

and by continuity 

ffoo(s (oo) ) = 1- 

Thanks to the continuity of solutions of stochastic differential equations with 
respect to both parameters and starting points (see e.g. Bahlali, Mezerdi 
and Ouknine pQ) we have, for any r > 0: 



lim P 

n— >+oo 



sup(X< n) (t) -X (oo \t)f 



t<T 



where X^°°^ is the solution of (S°°) starting from x^°°\ 

Therefore, for any A > 0, we infer from Lemma 15. II that if c T = Jq{P\{s) + 

p 2 N {s))ds for then 



lim P x 

n— >+oo 



(n) 



-A£(pf"(sf+p^(sf)ds 



P X 



-A£(p*(s)+p%(s))ds 



-Ac T 



Observe now that 
/o 



.7 



(p?( S )+p|K*))^ = £,H 



t s 



1 i 

2 _ fc 



(pf"(sf +p*"( S f)d S 



-En / 



Since k >2, choosing r = we see that there exists no such that toE% 
r for n > uq. Finally, 



1 _ i 

2 k \ 
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lim sup P x " 



Cf*°(pl(s)+p%(s))ds 



< lim sup P 
= 0. 



,(n) 



-CE^ j;{p^(sf+ P E N "(sf)ds 



Assume now that I < k. From the identity He„(x^) = 1 we cannot infer 
anymore that the sequence lives in a compact set. Nevertheless, by 
compactness, we can assume, by taking suitable subsequences, that 

• p- n ^(0), the Pi-th coordinate of x^ n \ converges to some for 1 < 
i < N. 

• 9i™^(0) ~~ converges to some 7$, for 2 < i < N. 

1_1 , n 

• 1 <7i (0) converges to a. 

Let be the process of coordinates {q^ l \t)—q^\fy,p^ l \t)). It is starting 
from x( n ) which lies in the hyperplane H = {x : q\ = 0}, and converges to 
some G K. The process iW is solution of the system 

(S E ") 

' dqi(t) = pi(t) dt 



dpi{t) = (-d qi H En - 
^d Pl (t) = (-d qi H En )dt 

with the Hamiltonian: 



E k 2 1 3 



- Pi ) dt + E* k 



TidBi(t) (i = l,N) 

(2 < i < N- 1). 



£*,(p,9)= E i^+^-^(^(ft-ffS n) (o)))+ E ^-^(ft+i-ft))- 



Ki<N 



E„ 



Ki<N-l 



Thanks again to the continuity of solutions of stochastic differential equa- 
tions with respect to both parameters and starting points we have, for any 
r > 0: 



lim P 



sup(A > ^(t) - X^°°\t)f 



t<T 







where X^ is the solution starting from x^ of the system (S°°). we can 
now finish the proof as in the case k = I. The only difference is that we have 
to prove that for a solution of (S 100 ) starting from x such that H QO (x) = 1 
and qi = 0, we have for any r > 0, 



(pj(s)+p 2 N (s))ds>0. 

Indeed, if we assume this integral to be 0, we obtain that on the time interval 
[0,t], the momenta are and the positions are constant : q(t) = c a vector 
solution of Vg(c) = if = ^ \p\ +g(c), so c is an infimum of the convex 
function g, and since c\ = 0, this entails c = which is a contradiction. □ 
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6. Proof of the main theorem 

We already know, from the preceding sections that the process is irreducible, 
strong Feller. To conclude the proof of TheorerdTTTJ we shall use Theorem 
8.9 of Rey-Bellet [9 J . There only remains to show that there exists to > 0, 
constants b n < +oo, < K n < 1, with lmin^+oo K n = 0, and compacts K n 
such that 

T t0 W(x) < K n W(x) + b n l Kn (x) 

If we choose compact sets of the type K n = {x : W(x) < a n } with a n — ► +oo, 
it is enough to show that 

hm sup — = 

{x-.w(x)>a n } yy{x) 



Assume that this is not the case, then we can find e > and a sequence x n 
with W{x n ) = e 0H ^ -» +oo and Tt ^f > e. 

Therefore, we shall finish the proof of Theorem ll.il once we have proved 

Lemma 6.1. Given < 9 < max(Ti,Tjv) , there exists a constant C > 
and a,(3>l conjugate exponents, ^ + 4 = 1, such that 



Zk*g*) < e ^(T 1+ T„)top,[ exp (-CJ\ P 1(S) + P 2 N (s))ds 



VP 



Indeed, by combining this Lemma with Lemma 15.21 we get that for a sub- 
sequence that we call x n again, 



e < e C9{Tl+TN)t ° limsu P P x " 
which is a contradiction. 



to 



exp -C / (pi(s)+p 2 N (s))ds 







Proof of Lemma \6.1\ The proof relies on the following clever trick (see e.g. 
Rey-Bellet and Thomas, proof of Theorem 3.10 [13]). Since L is the infini- 
tesimal generator of the diffusion X(t) = (q(t),p(t)) solution of (5), we have 
the decomposition 

H(X t ) = H(x) + f LH(X S ) ds + M H (t) 
Jo 

with M H (t) = J *cr'Vi7(X s ) dB s a continuous local martingale of quadratic 
variation 

r-t rt 

(M H ) = / ||cr*Vi/(X s )|| 2 ds = 2 / TH(X s )ds 
Jo Jo 

where T denotes the carre du champ operator associated to L: 

TH(x) = ^(L(H 2 )(x) - 2H(x)LH(x)) . 
Only the second order differentials play a role in T so 

TH(x) = \ Y, ^(d p ,(H*)-2Hd p .H) = \ Ti{d pi H? = \ £ T % p\ . 

i=l,N i=l,N i=l,N 



12 



PHILIPPE CARMONA 



Hence 
(6.1) 



LH{x) + a6TH{x) = -((Ti + T N ) - {p\ + p 2 N ) + a9{T W \ + T N p%)) 
< i(Ti +Tjv) - i(l - a0max(Ti,:ZV))(p? + 



Therefore 



W(x) 



px 
px 



J(H(X t )-H(x)) 



px 



9M{ I +8f*LH(X s )ds 



e eM»-a^{M») te a^(M») t+ ej*LH(X s )ds 



= P X [UV] <p^[u a ] 1 / a P x 
But U a is just an exponential martingale, 



V 13 



VP 



s a8M t H -±a 2 8 2 (M H ) t 



P x \U a ] =P a 



Hence, injecting the inequality (|6.ip . we obtain 



T tW{x) < i^ (Tl+Tjv) 
W(i) - 



x . . . 



x 



e 2 



(l-c«ema X (ri,T JV ))/ t (p?( S )+pl( S ))d S 



VP 



If we choose a close enough to 1 so that a0 max(Ii,Tjv) < 1, we obtain the 
desired inequality. □ 



7. The moment exchanging model 
The Hamiltonian is harmonic with no pinning (see [2]) 

l<i<N l<i<N-l 

Considering the interdistances = qi + i — we not only have two heat 
baths at both ends of the chain, but also a momentum exchanging noise. 
The stochastic dynamics are described by the stochastic differential system: 

(£) 

' dri(t) = (p i+1 -pi)dt 

dpi(t) = {n - — t2pi)tft - ^p 2 dB h2 {t) + y/T[dBo tl (t) 

1 + 7 / 

dpN{t) = (-rjv-i — Pn) dt + ^jp N _idB N _i iN (t) + ^T N dB N>N+1 {t) 

k dpi(t) = (rj - rj_i - TPi) + v / 7(pi-idB i _i )i (t) - p i+1 dB iii+ i(i)) (2 < i < JV - 1) 
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where {Bj i+i)o<i<jv are independent Brownian motions. The infinitesimal 
generator is now 

L = A + C R + |<S 

l<i<JV-l 

We have : X iji+1 (pf) = 2pip i+1 = -X iji+1 (pf +1 ). Therefore, X^+iH = 0, 
S(H) = and 

s(w) = w{e\x hl+1 Hf + exl +1 {H)) = o . 

i 

This implies that we can proceed through section 2 with no change. In 
particular: LW < \0(T X + T N )W if < 0m.ax.(T 1 ,T N ) < 1. 

Since we have more squared vector fields, we also have Hormander's hypoel- 
lipticity condition. We also dispose of more control, so the semigroup is 
irreducible (section 4). 

To use the scaling technique of section 5, we need to look closely at the 
semimartingale decomposition of H{Xt): 

H{X t ) = H(x) + Mf + f LH(X S ) ds 

Jo 

with 



dM t H = Pl {-^p 2 dB ly2 {t) + s/TxdBo^t)) 

+ ^ VlPi(Pi-l dB i-l,i ~ Pi+ldB iji+ i) 
2<i<N-l 

+ PN(VlPN-idB N - l!N + \/T^dB N , N+ i 

= pi y^cLBo,! (*) + PNdB N>N+1 (t) . 

Therefore, d(M H ^ t = (p\T\ + p 2 N Tjq)dt and the proofs of section 6 can pro- 
ceed (almost) unchanged. The only thing that changes is that the limiting 
system is no more deterministic but stochastic : only the heat reservoir 
noises disappear 

(S°°) 

' dri(t) = (p i+ i -pi)dt 
1 + 7 

dpi{t) = (n 2~pi) dt ~ -fyP2dB 1>2 (t) 

1 + 7 

dpN{t) = (-tn-1 y~ PN > dt + V^PN-idB N -i >N (t) 

x dpi(t) = (n - r;_i - 7^) dt + y/iipi-idBi-^iit) - p i+ idB iji+1 (t)) (2<i<N-l) 
We only need to replace Lemma f5.1l with the following non degeneracy result: 

Lemma 7.1. Let (q(t),p(t)) be a solution of (S°°), starting from x such 
that H OQ (x) = 1. Then, for any r > 0, almost surely, 

(pl(s)+ P %(s))ds>0. 
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Proof. Assume on the contrary that 

(pl(s) + p%(s))ds = a.s. 

o 

Then on / = [0, r], we have p\ = a.s., and thus 

i-t <-t 
ri (s)ds= / ^ P2 (s)dB h2 (s) , (Vi<r) 

'0 JO 

This states that a finite variation process coincides with a continuous mar- 
tingale. Hence, they both vanish : on I, r\ = p 2 = a.s. It should be 
now clear how to proceed by induction and to obtain rj = = Pi, which 
contradicts the fact that H (p, q) = 1 (here = H). □ 

8. The Lefevere-Schenkel chain 
Lefevere and Schenkel [7] consider a periodic lattice Hamiltonian 

1 N 

^ = 2 E^ + + ^ " qi -^ 

(we have pjy+k = Pk, QN+k = Qk and for sake of notations N is a multiple of 
4). They introduce the Fourier coordinates 

1 N 1 N 2 k 



i=i i=i 

which satisfy P_fc = (the complex conjugate) and = This 

implies in particular that Qo,Po, Qn/2, Pn/2 are rea l valued since for example 
Pn/2 = Pn/2-n = P-N/2 = Fourier coordinates the Hamiltonian 

reads: 

# = J E l^| 2 +u; 2 |Q fe | 2 



AT 



with |z| 2 = zz* the complex square modulus and uo\ = w 2 (/i 2 + 4sin 2 (^)). 
The momenta P k are coupled to heat baths at a fixed constant temperature 
T > 0, whereas the positions are coupled to heat baths at tempera- 
tures Du where to follow notations of Lefevere and Schenkel |7) we have 
D k = yT 2 a 2 (^). In short they are solutions to the system of stochastic 
differential equations 

, Lg , | dP k (t) = -(ulQ k + ^P k ) dt + VTdZ k (t) 

\dQ k {t) = P k dt - iD k dZ k {t) 

where {Z k {t))t>o are independent standard complex Brownian motions if 
k ^ 0,N/2 and ordinary real valued Brownian motions if k = 0,N/2. The 
boundary conditions imply that D k = P>~ k = Dk+N (and thus Dq = D N / 2 = 
0) and that Z_ k {t) = Z£(t). 

We shall work with real valued processes : P k = R k + iS k , Q k = V k + iW k , 
Z k = X k + iY k for k ^ 0, N/2. The stochastic differential system is for 
k ^ 0, N/2: 
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{LS A 



and for k 



(LS 2 



dR k (t) = -(u 2 k V k + -R k ) dt + VTdX k (t) 

dS k (t) = -{u 2 k W k + ~S fe ) dt + VTdY k (t) 

dV k (t) = R k dt + D k dY k (t) 
{ dW k (t) = S k dt - D k dX k (t) 

0, N/2 ( D k = and P k ,Q kl Z k are real valued): 
dR k (t) = -{u 2 k V k + ^R k )dt + VfdX k (t) , 



dV k (t) = R k dt. 



Therefore the infinitesimal generator can be written C 
Hamiltonian generator and 



A + C R with A the 



P=R k ,S k 



+ ^VTD k (dl 



14 



d R k W k ) 



Hormander's hypoellipticity condition is very easy to check and we have 

CH = AH + C R H = + NT + Y^ D 2 k ul - ^ |P fe | 2 < C , 

k k 

with C a constant, which we can choose C > 1, and with the function 
W = C + H that satisfies LW < W, we see that there is no explosion. 
The control problem is easy to solve since it splits into two elementary control 
and one in R 2 . Let us examine the control problem in 



problems, one in 

dx 
dt 



: 2| = Ax + Bu(t), with 



M 









( 


sfT 


\ 





1 

2 


o 


B = 







VT 


1 

















1 


J 




V- 


-D k 


o / 



A 



We let S be the smallest vector space containing the range of B and stable 
by A. S contains the vector f\ = (y/T, 0, 0, — D k ) T and Af± so it contains 
Afi + ^/i = (— oj 2 D k , 0, 0, VT) T and thus S contains the coordinate vectors 
(1, 0, 0, 0) T and (0, 0, 0, 1) T . Working similarly with f 2 = (0, VT, D k , 0) T we 
conclude that S is M 4 and by Kalman's criterion the accessibility sets are 
both M 4 . 

To use the scaling technique of section El we need to take a close look at 
the semimartingale decomposition of H(£t) with ^ = (P k (t), Q k (t), y + 1 < 
*<*). 



dH(£t) = dM t H + CH{i t ) dt 
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with M t a continuous local martingale: 

dM t H = Y,(VT - u 2 D k W k )dX k (t) + (VfS k + u? k D k V k ) dY k (t) , 

k 

with continuous quadratic variation d(M H ^ t = 2TH(^t) dt. Hence: 
TH(0 = \ ^(\/T - io 2 k D k W k f + (VfS k + u? k D k V k f 



2 

k 



<Y^T\P k \ 2 +u;iDl\Q 



2 

fcl 



Proceeding as in section [6J we see that 

LH(0 + 6aTH(£) < C + + a0r)|P fc | 2 + a^Z?||Q fc | 2 

Hence, if we fix a > 1, there exists constants #o> Ci> C2, C3 > such that 
for < < 6» , 

Lfr(o + ^r^(e) < d - ic 2 J] (|P fc | 2 - 9C 3 \Q k \ 2 ) . 

k 

This yields an analogue of Lemma 16,11 and the proof can proceed as in 
section [6l More precisely, we have the upper bound for < 9 < 9q , W = e 9H 



,1) T -§^<^p* 

W(x) 



1, 



V/3 



There is only one thing left to establish, the non degeneracy of the limiting 
system : the noise disappear, so that (R k , V k ) and (S k , W k ) are solutions of 



the same deterministic system, with to 2 = uii 




2 1 X , 

, 1 ~\0J v + -r)dt 

Lemma 8.1. Let x = (r,v) be a solution of with starting point xq 7^ 
(0,0). Then 

f T t 2 ( ds 
liminf J ° - ) ; > 0. 
j Q v^{s) ds 

Proof. This is a simple exercise in ordinary differential equations. The so- 
lution of the system is x(t) = e tA XQ. The eigenvalues of A are 



■| ± ii/cij 2 — ^ otherwise. 



Hence, t> = A + e x+t + A_e x ~t, r(t) = ^ and thus, in the first case, uj 2 < A 

/ T r 2 ( S )ds _ (A 2 , ifA^O; 



lim 



t^+oo J^v 2 (s)ds [A 2 otherwise. 



We leave the second case to 2 > ^ to the interested reader. □ 
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